Scattering matrix of elliptically polarized waves 



O 

5-1 

Oh- 
< 



^ ■ 

in 

o 

Of 



> 

oo 
m 

(N 

o 



X 



M. Martinez-Mares and E. Castano 
Departamento de Fisica, Universidad Autonoma Metropolitana-Iztapalapa, A. P. 55-534, 09340 Mexico D. F., Mexico 

We analyze the scattering of elliptically polarized plane waves normally incident at the planar 
interface between two different materials; we consider two cases: dielectric-dielectric and dielectric- 
conductor interfaces. The scattering matrix S in both cases is obtained using the boundary condi- 
tions and Poynting's theorem. In the dielectric-dielectric case we write S using two different basis, 
the usual xy and a rotated one. For the dielectric-conductor interface, the use of the rotated basis 
together with an energy balance argument leads us, in a natural way, to construct a unitary S matrix 
after recognizing the need to introduce two equivalent parasitic channels due to dissipation in the 
conductor, and the transmission coefficient into these parasitic channels measures the absorption 
strength. 

PACS numbers: 41.20.Jb, 42.25.Bs, 42.25.Fx, 42.25.Gy, 42.25. Ja 



I. INTRODUCTION 

The scattering of waves has attracted much attention 
even before being an object of scientific inquiry, since the 
observation of these phenomena is fascinating and even 
pleasurable. Take for example, the observation of water 
waves moving on the surface of a pond, where diffrac- 
tion and interference effects give origin very interesting 
patterns. Nowadays, scattering phenomena are use in 
general to study how a wave an a target are transformed 
by their mutual interaction, in such a way that we ex- 
tract very useful information about the structure of the 
target and the character of the wave itself. 1,2 

In wave phenomena, either classical or quantic, the dis- 
persion is mainly characterized, in an elegant and com- 
pact way, by the scattering S matrix, that describes the 
transformation of an incoming wave into an outgoing one 
due to the interactions with a particular target. 3 ' 4 This 
has been used to discover the inner structure of many dif- 
ferent objects ranging from macroscopic crystals, DNA 
molecules down to systems of atomic and nuclear sizes, 
and even smaller systems. 5 ' 6 

In electromagnetic phenomena scattering is produced 
by variations either in time or space of the dielectric 
function; 7 if this function is real then S is unitary since 
the energy flux is conserved. However, in the presence of 
a sink of energy the energy is not conserved as it happens 
in metals due to dissipation, in this case the dielectric 
function is complex and, therefore, S is a sub-unitary 
matrix. A large majority of the work on scattering has 
been done with unitary S matrices in problems where the 
flux is a conserved quantity, when dealing with problems 
where the flux is not conserved it then becomes natural 
to extend the S matrix to a unitary one by increasing its 
dimensions to include dissipation channels. 8-10 

In a previous publication, 10 we addressed the scatter- 
ing of linearly polarized plane waves with normal inci- 
dence at the planar interface of two media for two differ- 
ent cases, one with a dielectric-dielectric interface, and 
another with a dielectric-conductor interface. This was 
done using Poynting's theorem at the interface to define 



an S matrix; for the dissipative case, dielectric-conductor 
interface, the scattering matrix becomes sub-unitary, S, 
since the energy flux is not conserved. In this last case, 
we were able to write S as part of a unitary S matrix by 
the introduction of a single "parasitic channel" related to 
the energy dissipation in the conductor. However, that 
work only considered normal incidence of linearly polar- 
ized waves; in this work, we consider a more general case: 
elliptically polarized waves. 

In order to be self contained, in the next section we 
summarize the main ideas concerning Poynting's theo- 
rem. In Sect. Ill this theorem and the boundary condi- 
tions at the surface of two dielectrics to define a scattering 
matrix. Sect. IV is devoted to the dielectric-conductor in- 
terface where Poynting's theorem help us to extend the 
sub-unitary scattering matrix to a unitary one. 



II. ENERGY BALANCE EQUATION 

Assuming a harmonic time dependence of the electric 
E(r, t), magnetic H(r, t), and density current J(r, t) com- 
plex fields, the time-averaged Poynting's theorem is given 
by the real part of the equation 10 



where 



V-S P (r)=-ij(r).E*(r), 



S P (r) = -E(r) xH*(r) 



(1) 



(2) 



is the time-averaged Poynting's vector. The right hand 
side of Eq. (1) is the negative of the time-averaged work 
done by the fields, per unit volume per unit time, and 
represents the conversion of electromagnetic energy to 
thermal (or mechanical) energy. The integral of Eq. (1) 
over a volume V enclosed by a surface £ gives (only the 
real part is physically relevant) 



(3) 



2 



where. 

$ = <j> S P (r) • hda, (4) 
is the net flux $ of Sp(r) through S, and 

W=^/ J(r)-E*(r)dK, (5) 

is the time-averaged rate of work done by the fields if 
there are dissipative processes in the system. 

A. Energy flux conservation 

When there is no dissipation, J • E* =0, and therefore 
W is zero, and the energy flux is conserved: 



$ = 0. 



(6) 



This means that the net flux crossing into the system 
equals the one leaving it. 



III. DIELECTRIC-DIELECTRIC INTERFACE 
A. Linear xy polarization basis 

Lets take the xy-plane as the surface that separates 
two dielectrics with indices of refraction n and n' ', as 
shown in Fig. 1. We consider a normally incident plane 
wave which is the superposition of two linearly polarized 
waves, one in x-direction and the other one in y, what 
we call x — y basis. 

Therefore, the spatial part of the electric field for z < 
is given by 



E(z) = {E ax ± + E ay y) e lkz + {E hx ± + E by y) e 



-ikz 



(7) 



where k — nu> /c; the subindex a denotes incoming waves 
and b outgoing ones. Since we are using an elliptically 
polarized wave we have that 



E ax = aE a and E ay = (3 E a e* 



(8) 



where <p a is the phase difference between x and y com- 
ponents, being a and (i real numbers whose squares add 
up to one, 



1; 



(9) 



E a is the magnitude of the electric field of the incident 
plane wave for z < 0, 

E a =E ax x + E ay y; (10) 

Similarly, on the right hand side, z > 0, we have that 

E'(z) = (E' bx ± + E' by y) e* k '* + « x x + E' ay y) e^ k ' z 

(11) 




FIG. 1. Elliptically polarized waves normally incident at the 
planar separation between two dielectrics with indices of re- 
fraction n and n' , respectively. 



where k' = n'uj/c and 



E 'ax 



a'E' a and E' = 0' E' a e^ 



(12) 



where c/)' a is the phase difference and E' a is the magnitude 
of 



El 



E' ax ± + E'y, 



(13) 



and a' and /?' satisfy an equation equivalent to (9). 
Equations (8) and (12) can be written in a matrix form 



E 

E ax 



where 
D n = 



= D n 



a 
a' 



E a 

EL 



and 



and Dg = 



E 

J->ay 
E 'ay 







= D< 





/3'e i4 



E a 
EL 



(14) 
(15) 



that, in correspondence with Eq. (9) must satisfy the 
following condition 



Di + DlD p = Di + D D} = l 2 . (16) 
What do next is to find the electric fields of the out- 



it 



going waves, 



E fc 
El 



E bx x + E by y, 



E' bx i + E' by y, 

in terms of the incoming ones. 

Imposing the boundary conditions at the interface, z 
0, we obtain Fresnel equations in a matrix form 1 



(17) 
(18) 



/ E bx \ 
' E! > 



bx 



S F 


2 


o 2 


S F 



( E ax ^ 

E'ax 



Eay 
\KyJ 



(19) 



where 2 is a 2 x 2 null matrix and Sf is a 2 x 2 matrix 
given by 



-r F t' F 
t F -r' 



(20) 
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where, rp (r' F ) and tp (t' F ) are the reflection and trans- 
mission Fresnel coefficients for electric field waves incom- 
ing from the left (right), whose explicit form, in terms of 
the indices of refraction, are 7 



rp 



t F = 



n' + n 

2n 
n' + n 



,r F = -r F , 

_ 2n> 
' tp ~ n' + n 



(21) 
(22) 



Following Ref. 10 we introduce renormalized electric 
field amplitudes given by 



a, b; m 



£i m = VnE lm , £[ m = VriE' lm (I 
Therefore, Eq. (19) can now be rewritten as 



where 



/ £bx \ 




( &ax \ 


c bx 




'-'ax 


£by 


Sxy 


£ay 






\£'ayj 




_ ( S 2 


o 2 \ 


Sxy 


\ o 2 


s 2 )' 



x,y). 
(23) 



(24) 



(25) 



with 




(26) 



being the 2x2 scattering matrix for normal incidence of 
a single linearly polarized plane wave (it was named S in 
Rcf. 10). 

The scattering matrix S xy defined in Eq. (24), is a 
4x4 unitary and symmetric matrix because S 2 is itself 
a 2 x 2 unitary and symmetric matrix. 10 Therefore, flux 
conservation is fulfilled as well as time inversion invari- 
ance; S xy is block diagonal since it is written in the basis 
of linear polarization, where the x and y components, 
called channels in the nomenclature of nuclear physics, 
are decoupled from each other. However, at this stage 
is not completely clear how to recover the simpler result 
of a single linearly polarized wave; we remedy this by a 
change of basis as shown in the next subsection. 



Rotated basis 



From Eqs. (14) and (24) we have that 
/ £bx \ 



°bx 


_ ( S 2 D a 


2 


£by 


V o 2 


S 2 D 



(£a\ 

£a 
\£'aj 



(27) 



Now, if a 7r/4-rotation is applied to the basis used up to 
now, by means of 13 



Ro 



J_ / t 2 1 2 
V2 \ -1 2 1 2 



(28) 



equation (27) becomes 



/ £bi \ 




( V2£ a \ 


c bl 


= S' 


V2£'a 


£bi 











V / 



(29) 



where we have defined 



£bl 
°bl 



<--b2 
£'b1 



1 / £by 

1 / £by 
V2\£'by 



_j_ / £bx 

V2\£L 

_j_ / £bx 

V2\£L 



and 



S' = Rq 



S 2 D a 


2 


2 


S 2 D I3 



R 



■0 1 



(30) 
(31) 

(32) 



is a new scattering matrix. It is important to realize that 
S' does not have the familiar form of a scattering matrix, 
as can be see by (29). 

In a standard scattering matrix the 11 and 22 blocks 
relate the incoming amplitudes on one side to the outgo- 
ing ones in the same side, and the 12 and 21 blocks re- 
late the incoming amplitudes on one side to the outgoing 
ones on the other side. Since we want to keep the famil- 
iar and very useful interpretation of a scattering matrix 
we transform (29) by means of the following orthogonal 
transformation matrix 



O 



( I o\ 

10 
10 

\ 1 j 



(33) 



that rearranges the incoming vector in a way that can be 
more easily interpreted in scattering theory. Therefore, 
we obtain that 



(34) 



/ £bi\ 




( V2£a\ 


£b2 


= s 





£U 


V2£'a 






K J 



where 



S = OR 



S 2 D a 


2 


2 


S 2 D P 



R T O T , 



(35) 
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is a scattering matrix whose dimension is twice the num- 
ber of channels. This scattering matrix can now be di- 
rectly identified as: 



S = 



r a 

t r' 



(36) 



where r (r') and t (?) are immediately identified as the 
corresponding 2x2 reflection and transmission matrices, 
respectively, for incidence from the left (right); they are 
explicitly given by 



n — n' 2\/nn' „ 

n + n' 
2 V / W 



n + n 
n — n 



n + n' 



F' t 



n + n' 



F' 



where 



1 

F =2 



F' = 



pe^" + a fie**" - a \ 
^ I3e 1 ^ - a fie^o + a J 

1 / P'e^'- + a' P'e 1 ^ - a' 



p'e 1 ^ - a' p'e^ + a' J 
Even though S satisfies flux conservation, 

s f s = u, 

docs not possesses a time reversal invariance, 



(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



where the superscript T denotes transposition. Also, 
since there is no specular symmetry we have that r' =/= r. 

We ask the reader to interpret why S xy is symmetric 
while S is not. (Suggestion: see the discussion leading to 
the derivation of Stoke's relations. n < 12 ) 

The reflection and transmission coefficients are now 
given by 



v ' {n + n') 1 



tr (i?) 



Ann' 
(n + n 



l\2- 



(43) 
(44) 



where we used that tr(FF^) — tr(F'F'^) = 1, where we 
note that due to flux conservation 



R + T = l. 



(45) 



Given all this, we can write the electric fields, Eqs. (30) 
and (31), as 



E'(z) = ^Se ie - 4fc ' z 



ei 



b =le 2 ) e~ ik * 
n 



(46) 



^ke, + ^Le 2 ) e ik '\ (47) 



where we have introduced a new basis set of vectors 

ei = -^=(y + x) and e 2 = -^=(y-x), (48) 

which plays a role equivalent role to the basis used in 
Ref. 7 to discuss circular polarization. 

We left to the reader the analysis necessary to verify 
that S reduces to the linear polarization case studied in 
Ref. 10, as well as to study the circular polarization case. 



IV. DIELECTRIC-CONDUCTOR INTERFACE 

We now study a different case. Lets assume that the 
media on the right side of the interface is a conductor 
with an electric conductivity a, such that its refractive 
index is complex: n' — > n' + in', where n' and n' are the 
optical constants; 14 the corresponding wave number is 
also complex. In the treatment of last section we replace 



where 



k' = n'uj/c and k' = t]'uj/c. 



(49) 



On the dielectric side, from Eqs. (46) and (47) with 
£' a = we have that 

E (z) = ^e ie fc + ( %e! + %e 2 ") e~ ik * (50) 



E'(z) = 



f 



61 



f» 1 _|_ 

\pn! + in' \/n' + in' 



f 

c b2 



e- 2 i e K 'v fe ' z . 



(51) 



By definition, the scattering matrix relates the outgo- 
ing to the incoming plane wave amplitudes; therefore 




where 



S 



n — n 



V2£ a 



in 



n + n' + in' 



(52) 



(53) 



Note that now S describes the reflection back to the same 
side where the elliptically polarized planes waves are ar- 
riving in a normal direction to the interface; therefore, it 
is a 2 x 2 matrix with the following structure 



S = 




(54) 



where now r (f 1 ) is the reflection amplitude when inci- 
dence is on the channel ei (e 2 ), and t (?) is the transmis- 
sion amplitude from channel ei (e 2 ) to channel e 2 (ei). 
Here, however, ei is the only channel where incidence is 
possible; note also that ? = i such that 



S = S 1 , 



(55) 



■5 




FIG. 2. On the right hand side, z > 0, there is a conductor 
with conductivity o, and for z < we have a dielectric mate- 
rial. The arrows represent Poynting's vectors corresponding 
to the electric fields on the dielectric. 



since there is time reversal invariance. Since channels ei 
and e2 are equivalent to each other, there is reflection 
symmetry thing for which f' — f. It is very important to 
realize that S is not unitary, and that according to Eq. 
(53), can be parametrized as 



S = y/Re l6 F, 
where R is the reflection coefficient given by 



R 



(n-n') 2 +n' 2 
(n + n'f + rj' 2 ' 



(56) 



(57) 



and 9 is a measure of the phase shift between incoming 
and reflected waves, and satisfies 



tan 6 = 



-2ni]' 



n' 2 - rf ' 



(58) 



Lets see how an energy balance consideration impose 
restrictions on S. Consider a closed surface £ consisting 
of a semi-infinite cylinder of cross section A, as shown 
in Fig. 2; in this case, Poynting's vectors are directed in 
the z-direction and the energy flux takes place through A 
and not through the lateral surface. The energy balance, 
Eq. (3), is written as 



$bi + $b2 - $ a = -W, 
where, based on (2) and (4), 10 



= 



1 



2/x c 
1 

2^oc 



\£ bj \ 2 A, j = l,2, 



V2£a 



A, 



(59) 

(60) 
(61) 



From Eq. (5), the rate of work made by the field is then 
given by 

W=\T ° (\£' bl f + \£' b2 f) e- 2 «'Mz, 
2 Jo \Jri 2 + rf 2 v 1 

(62) 



where we use that 3(z) — (tE'(z); integrating we have 
that 

W = - , a fe| a + |£k| 2 )- (63) 

Using Eqs. (60), (61) and (63), Eq. (59) gives 

(|£ hl | 2 + |£ b2 | 2 )-2|£ a | 2 = ~ „ T, C ,M ,? + \^\% 

(64) 



2n\/n' 2 + r]' 2 
which on a matrix form is written as 



(V2£* o)(5 f S-l2) 



V2£a 





(7fl C 



2ny/n' 2 + 



( cl* CI* \ I 

—pt \ & 61 £ M I ci 



I 1 

^b2 



V2£a 



where [see Eqs. (37)] 



t _ 2^/n{n' + irf) F 
n + n' + it]' 

Therefore, Eq. (65) is equivalent to 
^5 + 4^ = 12, 

where 



in 



j 2 K y/n' 2 + ri' 2 
which can also be written as 



t. 



where now 



and <fi is given by 



Ann 1 



(n + n' f + rf ' 



tan0 



rf 



n + n' 



(65) 



Now, using Eq. (34) with £' a — and n' — > n' + in' we 
can write that 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



Calculating the trace of Eq. (68) we obtain an expres- 
sion that in short is equivalent to energy flux conservation 



R + T v = l, 



(73) 



where R, the reflection coefficient is the fractional 
amount of energy that is turned back into the dielectric, 
while T p tell us the amount of energy dissipated or "lost" 
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in the conductor to non propagating modes or, what we 
call parasitic channels. 9 ' 10 

Hence, Eq. (68) can be seen as the unitarity condition 
for a 4 x 4 scattering matrix S that satisfies flux conser- 
vation, and is given by 

S = ( S tp J , (74) 

where the unitarity of S implies that 

S pp = -\~Sh p = -^T^%e^-^F. (75) 
tp 

The lack of unitarity of S is then given by following ex- 
pression 

P = I 2 - &S = th p , (76) 



whose trace quantifies ohmic losses through an absorption 
strength parameter 

7 = trP = T p , (77) 

where it must be emphasized that an ohmic loss takes 
place through parasitic channels in the language em- 
ployed in scattering theory. 

Therefore, the S'-matrix given by Eq. (74) describes 
how the energy losses in the conductor can be interpreted 
as energy lost to parasitic channels characterized by the 
t p matrix, where the coupling between the two media 
is described by the material constants of our model; in 
particular, in this work, for the elliptic polarization case. 
Even though we have two parasitic channels the total ab- 
sorption is measured by a single parameter 7 = T p , that 
gives the same result previously obtained in the linear po- 
larization case, something expected since each polariza- 
tion mode contributes just a fraction of the total energy 
flux. 
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